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Abstract 

In this paper, we study wave properties of isothermal plasma for 
the Schwarzschild de-Sitter black hole in a Veselago medium. We 
use ADM 3 + 1 formalism to formulate general relativistic magneto- 
hydrodynamical (GRMHD) equations for the Schwarzschild de-Sitter 
spacetime in Rindler coordinates. Further, Fourier analysis of the lin- 
early perturbed GRMHD equations for the rotating (non-magnetized 
and magnetized) background is taken whose determinant leads to a 
dispersion relation. We investigate wave properties by using graph- 
ical representation of the wave vector, the refractive index, change 
in refractive index, phase and group velocities. Also, the modes of 
wave dispersion are explored. The results indicate the existence of the 
Veselago medium. 
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1 Introduction 



Our solar system is filled with a wide range of celestial objects. Black hole 
is one of such objects, having so strong gravitational pull that no nearby 
matter or radiation, not even light can escape from its gravitational field. 
Astronomers are curious to extract real life examples of black hole. The 
presence of matter in the form of white dwarfs and neutron stars suggests 
the existence of stellar mass black holes. The accumulated evidence for the 
black hole existence is now very captivating. It is believed that collapse of 
a massive star under its own gravity leads to the formation of black hole 
(Das 2004). Plasmas are abundant in nature, almost found everywhere in an 
interstellar medium. It is a distinct state of matter with free electric charge 
carriers which behave collectively and respond strongly to electromagnetic 
fields (Raine and Thomas 2005). Black hole (in its surroundings) attracts 
plasma towards the event horizon due to its strong gravitational pull. The 
plasma forms an accretion disk due to interaction of plasma field with black 
hole gravity. 

The theory of general relativistic magnetohydrodynamics (GRMHD) is 
the most reliable discipline to examine the dynamics of magnetized plasma 
and effects of black hole gravity. The de-Sitter spacetime is a vacuum solu- 
tion of the Einstein field equations including a positive cosmological constant 
(Rindler 2001). The Schwarzschild de-Sitter (SdS) metric describes a black 
hole expressing a patch of the de-Sitter spacetime. Since the SdS black hole 
is non-rotating, so plasma in magnetosphere moves only along the radial di- 
rection. According to the recent cosmological and astrophysical observations, 
our universe is accelerating rather than decelerating and inclusion of positive 
cosmological constant reveals the expanding universe (Reiss et al.1998; Bah- 
call et al. 1999; Perlmutter et al. 1997). That is why our universe approaches 
to de-Sitter universe in future. This motivates the study of plasma waves in 
de-Sitter spacetime. 

Petterson (1974) investigated the strong gravitational field close to the 
surface of compact objects for the Schwarzschild black hole. Narayan (2005) 
suggested that compact objects having mass three times the solar mass can 
be identified as black hole candidates. Plasma present in magnetosphere 
is perturbed by gravity of black hole. Zerilli (1970a, 1970b, 1970c) used 
linear perturbation to explore gravitational field of a particle falling in the 
Schwarzschild black hole. Price (1972a, 1972b) discussed dynamics of ap- 
proximately spherical star by using non-spherical perturbations. Regge and 
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Wheeler (1957) also used non-spherical perturbation to investigate the sta- 
bility of Schwarzschild singularity. Gleiser et al. (1972) explored the stability 
of black holes by considering second order perturbations. 

Arnowitt, Deser and Misner (ADM) (1962) proposed 3 + 1 formalism for 
an easy approach to General Relativity (GR) by separating metric field into 
two parts (space and time) to characterize the coordinate system. Smarr and 
York (1978) used this formulation to explore spacetime kinematics numeri- 
cally. Israel (1967, 1968) discussed event horizons in static vacuum and static 
electro-vacuum spacetimes. Thorne and Macdonald (1982a, 1982b) explained 
how 3 + 1 split is appropriate approach for black hole theory. Macdonald and 
Suen (1985) developed a self-consistent formalism to treat electromagnetic 
and gravitational fields near black hole horizon. Sakai and Kawata (1980) 
analyzed wave propagation in ultra-relativistic plasma, parallel to a constant 
magnetic field in a frame of two fluid model. Holcomb (1990) and Dettmann 
et al. (1993) constructed electrodynamical equations for the universe mod- 
els. Holcomb and Tajima (1989) formulated linearized theory for relativistic 
plasma and found results for matter fluctuations in the early universe. 

Rezolla et al. (2003) explored dynamics of thick disks around SdS black 
hole by considering the effects of cosmo logical constant. Font and Daigne 
(2002) studied stability of thick accretion disks around black holes. Myung 
(2001) developed entropy bounds for SdS black hole. Suneeta (2003) con- 
sidered quasinormal modes for scalar field perturbations of SdS black hole. 
Setare (2005) obtained area and entropy spectrum near extremal SdS black 
hole horizon. Zhang (1989a) modified the stationary symmetric GRMHD 
black hole configuration theory. The same author (Zhang 1989b) explored 
the modes of perturbation in rotating black hole. Buzzi et al. (1995a, 1995b) 
determined the properties of waves propagating in two fluid plasma for the 
Schwarzschild black hole. Ali and Rahman (2009) explained transverse wave 
propagation in two fluid plasma around SdS black hole. Sharif and his collab- 
orators (Sharif and Sheikh 2007a, 2007b, 2007c, 2008a, 2008b, 2008c 2009a, 
2009b; Sharif and Mustafa 2008; Sharif and Rafique 2010) have explored 
wave properties of cold, isothermal and hot plasmas with Schwarzschild as 
well as Kerr spacetimes in the usual medium. 

The medium with both negative permeability and permittivity has the 
unusual electromagnetic properties named as Veselago medium or negative 
index medium (NIM), after a Russian physicist Veselago (1968). It is also 
called as double negative medium (DNM) or negative phase velocity medium 
(NPV). Valanju et al. (2002) presented treatment for refraction of electro- 
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magnetic waves in a NIM. Ross et al. (2006) concluded that propensity 
of a rotating black hole is enhanced in the presence of charge to support 
wave propagation with negative phase velocity in its ergosphere. Ziolkowski 
and Heyman (2001) studied wave propagation analytically and numerically 
in NIM. Nagar et al. (2004) reported results from numerical simulations of 
gravitational radiations emitted due to matter accretion from non-rotating 
compact objects. In recent papers, Sharif and Mukhtar (2011a, 2011b) have 
discussed wave properties with non-rotating as well as rotating background 
plasmas (isothermal and hot) in this unusual medium. 

This paper deals with wave properties of isothermal plasma around SdS 
black hole in a Veselago medium. We consider 3 + 1 GRMHD equations and 
determine a dispersion relation by Fourier analysis for both magnetized and 
non-magnetized backgrounds. The results are discussed by three dimensional 
plot of wave vector, refractive index and change in refractive index. The pa- 
per is organized as follows: In section 2, linearly perturbed 3 + 1 GRMHD 
equations for isothermal plasma and their Fourier analysis is developed. Sec- 
tions 3 and 4 provide reduced form of the GRMHD equations for rotating 
(non-magnetized and magnetized respectively) plasmas. We summarize our 
results in the last section. 

2 GRMHD Equations in a Veselogo Medium 
With Isothermal Plasma Assumption 

The general line element in ADM 3 + 1 formalism is given as follows (Zhang 
1989b) 

ds 2 = -a 2 dt 2 + rnjidx* + ^dt)(dx j + frdt). (2.1) 

A natural observer associated with this spacetime is known as fiducial ob- 
server (FIDO), a denotes lapse function (ratio of FIDO proper time to uni- 
versal time i.e., ^), (3 % is three-dimensional shift vector (which determines 
change in spatial coordinates) and 77^ = 1,2,3) are the components of 
three-dimensional hypersurfaces. The SdS spacetime in Rindler coordinates 
is given by (Ali and Rehman 2009) 

ds 2 = -a 2 (z)dt 2 + dx 2 + dy 2 + dz 2 , (2.2) 

where the directions z, y and x are analogous to the Schwarzschild coordi- 
nates r, (f) and 6 respectively. Since SdS black hole is non-rotating, the shift 
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vector vanishes. On comparing Eqs. (l2.1| ) and (\2.2\i , we have 

a = a(z), f3 = 0, 77^ = 1 (i = j). (2.3) 

The 3 + 1 GRMHD equations for the line element (12.21) in a Veselago 
medium are given by Eqs. (1A5l )-( lA9l) in Appendix. The equation of state for 
isothermal plasma is (Zhang 1989a) 

p = = constant, (2.4) 

Po 

here po, p, p and p denote rest mass density, moving mass density, pressure 
and specific enthalpy respectively. The specific enthalpy is constant while 
pressure p 7^ for the isothermal plasma. This equation shows that there 
is no energy exchange between plasma and magnetic field of fluid. The 
corresponding 3 + 1 GRMHD equations ( (IA5j) - flA9j) ) for isothermal plasma 
surrounding the SdS black hole become 

_ = _Vx (qVxB), (2.5) 
V.B = 0, (2.6) 
fZe + P). + ( p + p) 7 2 V.^ + (aV.V)(p + p) + (p + p) 7 V.(aV.V) V 
+ (p + p)V.(aV) = 0, (2.7) 

(p + Ph 2 + ?T) % + (p + Ph 4 ^- - ^ W 



+V.V) V* - ( — 5 V - —BiBj V; k V k 



B 1 

— &a 

4tt 3 4vr 

+ ( p + p) 7 2 a . + = J_ (v x B ),V.(V x B) - -^(aB)5 

47r 87rcr 

+ i -L(a£ i ) J i^ - _L[B x {V x (V x (oV x B))}] if (2.8) 

{— + V.V)(p + p)f - ~ + 2(p + p)f(V.a) + (p + p) 
a or a at 

7 2 (V.V) (V x B).(V x ^) - — - (V x B).(B x ^-) 

+— (VxB).(VxaB)=0. (2.9) 
Ana 
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In rotating background, we assume that plasma flow is in two dimensions, 
i.e., in xz-plane. Therefore FIDO's measured velocity V and magnetic field 
B turn out to be 

V = V(z)e x + u(z)e z , B = B[X(z)e x + e z ], (2.10) 

here B is an arbitrary constant. The relation between the quantities A, u 
and V is given by (Sharif and Sheikh 2007a, 2007b, 2008a, 2008b, 2008c) 

V F 

V = + \u, (2.11) 

a 

where V F is an integration constant. The Lorentz factor, 7 = , 1 be- 
comes 

When the plasma flow is perturbed due to black hole gravity, we use 
linear perturbation. The flow variables (mass density p, pressure p, velocity 
V and magnetic field B) take the form 

p = p° + Sp = p° + pp, p = p° + Sp = p° + pp, 

V = V° + SV = V° + v, B = B° + 5B = B° + Bb, (2.13) 

where unperturbed quantities are denoted by p°, p, V° and B°, the linearly 
perturbed quantities are represented by Sp, Sp, 5V and SB. We introduce the 
following dimensionless quantities p, p, v x , v z , b x and b z for the perturbed 
quantities 

p = p(t,z), p = p(t,z), v = SV = v x (t, z)e x + v z (t, z)e z , 

b = = b x (t, z)e x + b z (t, z)e z . (2.14) 

When we insert these linear perturbations in the perfect GRMHD equations 
(Eqs.( f2T5i) -f l2.9p ) along with Eq. fl2.14p . the component form of these equations 
will be (Sharif and Mukthar 2011a, 2011b) 

~~^T ~ uh *,z = ( ub *= ~ Vb * ~ v x+ ^z) V In a 
a at 

-K* - Xv ZiZ - X'v z + V'b z + Vb Z;Z - u'b x ), (2.15) 

1 db z , 

-7^ = °> 2 - 16 
a at 

b z , z = 0, (2.17) 
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dp dp . ,,,3^ dv z . 

P ~di + P ~di + ( p + p W ( ~dt + U ~dt' + aupp ' z + aupp > z 

+a(p + p){'~i 2 uVv X)Z + (1 + ^ 2 u 2 )v Z)Z ] - -(p-p)(au^p) tZ 

7 

+a(p + p)7 2 w{(l + 2~f 2 V 2 )V + 2~f 2 uVu'}v x - a(p + p) 

x{(l - 2 7 V)(1 + 7 V)^} - 2 1 \ 2 VV'}v z = 0, (2.18) 

(p + ph 2 (i + 7 V 2 ) + 1% + ((p + P )Av - A /] 



4ir I a 9t I 47r 

1 <% 2 



a <9t 



+ | (p + p) 7 2 (l + 7V 2 ) + £- J m^, 2 + { (p + p) 7 4 ^ 

^ ~ + w2)6 "' 2 ~ 4^ + U ^ + W } ^ 
+7 2 w(pp + pp) { (1 + 7 2 y 2 )y' + ^uVu'} + [(p + p)i A u{(l 

B 2 va' 

+A 1 2 V 2 )uu' + 4W(i + 7 2 v 2 )} + — —]v x + [(p + p) 7 2 {(l 

47TQ: 

+2 7 V)(1 + 2 7 V 2 )\/' - 7W + 2 7 2 (1 + 2 7 2 m>F«'} 

-^ (Aa) ^ = ' (2 ' 19) 

A 2 B 2 1 ldu, r, , a „ AS 2 



(p + p) 7 2 (l + 7 V) + — - -— ^ + (p + p) 7 4 ^ 



An ) a dt { ' An 

x ^ ^ + I (p + P)l 2 (1 + 7 V) + } + { (P + Ph 4 u 

xV- ^-}uv XiZ + ^-(1 + m 2 )6 Xi2 + -^-{c/A - (a\)' + uX 
An An Ana 

x(ua' + u'a)}b x + (pp + pp) 7 2 {a 2 + uu'(l + 7V) + j 2 u 2 VV'} 
+ [(P + ph 4 {u 2 V'(l + 47 2 ^ 2 ) + 2F(a 2 + uu'(l + 2 7 V))} - AS 2 

X T^K + [(p + P)7 2 {«'(1 + 7V)(1 + 47 2 ^ 2 ) + 2w 7 2 (a 2 + (1 
Ana 

+2 1 2 u 2 )VV')} + ^^(aA)'K + (p'p + Pp') = 0, (2.20) 
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1 dp 1 dp 

~^Plu + ~^ 2 Plu + + P'H + U l 2 (PP,z + PP,z + p'p 

a at a at 

+p'p) P^r + 2j 2 u(pp + pp)a z + j 2 u'(pp + pp) + 2(p 

a or 

+p)7 4 (uV r/ + 2uVa z + mV)?;, + 2(p + p) 1 2 (2 1 2 uu + a 27 4 
+2 7 Va> 2 + 2(p + p) 7 4 wl^ + (p + p) 7 2 (l + 2 7 2 w> 2 , 2 

B2 [(V 2 + u 2 )X^ + (y 2 + , 2 )^ - F(AF + u) dK 



4vra L <9r 7 dt dt 

db B 2 

-u(XV + u)-^] + —[(V - \u)v Xtt + X(uX - V)v g7t ] 
B 2 

-—(XX'v z - \'v x - X'Vb z + X'ub x - Vb xz + uXb x ^ z ) = 0. (2.21) 

The following harmonic spacetime dependence of perturbation is assumed 
for the Fourier analysis, 



p{t, z) 


= c\e 


-L(L)t- 


-kz) 

i 


p(t, z) 


= c 2 e- i{wt - kz \ 


v z (t,z) 


= c 3 e" 


-i(uit- 


-kz) 

; 


v x (t,z) 


i = c 4 e- L{ult - kz) 


b z (t,z) 


= c 5 e" 


-i(uit- 


-kz) 


b x (t,z) 





(2.22) 

here u and k represent the angular frequency and ^-component of the wave 
vector (0,0, k), respectively. The wave vector can be used to determine re- 
fractive index and the properties of plasma near the event horizon. 

The ratio of speed of light when it travels from one medium to another is 
said to be refractive index. Frequency dependence effects in wave propagation 
refers to dispersion. This describes relations between wave properties like 
wave length, angular frequency, refractive index etc. (Das 2004). Dispersion 
is said to be normal if change in the refractive index with respect to angular 
frequency is positive, otherwise anomalous. Using Eq. fl2.22p in Eqs.f l2.15j) - 
f)2.2ip . we get their Fourier analyzed form 

Ci{a' + ika) — C3 {(aX)' + ikaX } — c^(aV)' + c§{(au)' + loj 
+ikua} = 0, (2.23) 

cs(— ) = 0, (2.24) 
a 

c 5 tk = 0, (2.25) 



8 



Ci{(— lui + ikau)p — pj 2 au(VV' + W) — a'wp — au'p — aup'} 
+c 2 {(— uuj + ikau)p + a'-up + cm'p + crop' + p^ 2 au(VV' + W)} 
+c 3 (p + p)[-io;7 2 M + tka(l + 7V) - a{(l - 2>y 2 u 2 )(l + 7V) 
u' 

x — - 2-f A u 2 VV'}] + c 4 (p + p)[j 2 V(-loj + ikau) + aj 2 u{(l 

+2 1 2 V 2 )V + 2 1 2 uVu'}} = 0, (2.26) 
c lP1 2 u{(l + 7 2 y 2 )y' + -y 2 uVu'} + c 2 p 7 2 M{(l + 7 2 y 2 )y' 

A -B 2 

+7 2 wW} + c 3 [-{(p + p)7 4 «F — }— + {(p + p)7 4 wV 

47r a 

+ (p + p) 7 2 {(l + 2 7 V)(1 + 2 1 2 V 2 ) - ^V^V 

An 

B 2 u 

+2 7 4 (p + p)uVu'{\ + 2 7 V) - =-%A)'] + c 4 [-{(p + p) 7 2 

Ana 

(1 + TV 2 ) + f }- + {(p + p) 7 2 (l + 7^ 2 ) + f}tku 
An a An 

+{p + ph 4 u{(l + 4 7 2 V 2 )W + 4W'(1 + 7V 2 )} + 

-c 6 ^{(l + w 2 )^ + (1 + u 2 )- + W} = 0, (2.27) 
An a 

cip7 2 {a 2 + uu'{\ + 7 2 m 2 ) + 7 2 «V7'} + c 2 [p7 2 {a 2 + uu 

(1 + 7 V) + 7 2 uW} +p' + ikp] + c 3 [-{(p + p) 7 2 (l + 7 V) 

A 2 £> 2 to; A 2 S 2 

+ -—} — + {(p + P)7 2 (l + 7 2 ^ 2 ) + -r-}^u + {(p + p) 7 2 
47r a An 

x{u'(l + 7 2 m 2 )(1 + 4 7 V) + 2m 7 2 {o z + (1 + 2 1 2 u 2 )}VV'} 

+^(a\)'} + c 4 [-{(p + p) 7 4 ^ - + {(p + p) 7 4 ^ 

47ra 47r a 

+ {(p + p)7 4 {«V(l + 4 7 2 \/ 2 ) + 2V{(1 + 2 7 V)W 

47T 

XB 2 a'u B 2 

+a z }] + Cq[- {— (aA)' + c/A — wA(mq; / + u'a)} 

Ana Ana 

\B 2 

+-— {l + u 2 )ik} = 0, (2.28) 

47T 
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ci{(— —J 2 + + 2u~f 2 a z + ~f 2 u')p + up'-f 2 } + c 2 {( — (1 - 7 2 ) 

a a 

+tku~f 2 + 27 2 wa z + ~i 2 u')p + W7V} + c 3 7 2 {(p' + p') + 2 

A.B 2 

x (2 7 2 W + a z + 2^u 2 a z ){p + p) + (1 + 2 7 V)(p + p)tk + 

47TQ; 

x (Aw - + a\'} + c 4 [2(p + p)^{{uV + 2uVa z + u'V) + uVik) 

d2 — R 2 

+- — (V - u\)lu - aX'} + c 6 [- {{V 2 + u 2 )X + V(XV + u)iu} 

vkol Aita 

-aX'u + ika{V - uX)) = 0. (2.29) 



3 Rotating Non-Magnetized Flow 

For the rotating non-magnetized background of plasma flow, we substitute 
B = = A and C5 = = cq in the Fourier analyzed perturbed GRMHD 
equations ( (T2T26l) -( l2T29l) ) (Sharif and Mukhtar 2011a, 2011b). 

3.1 Numerical Solutions 

For the rotating non-magnetized plasma, we use the following assumptions 
to find out the numerical solutions 

1. Specific enthalpy: p = 1, 

2. Time lapse: a = 

3. Velocity components: u = V, x and ^-components of velocity yield 



u = V 



1 



4. Stiff fluid: p = p= 

5. Lorentz factor: 7 



7751 



Vl-u 2 -V 2 2 ' 

where is the SdS event horizon greater than that of the Schwarzschild 

event horizon and r h « 2M (l + ^f- + - C0.2948A;m, 1 ^ £ ^ 1.5 for a 

black hole mass M ~ IMq. The value of £ corresponding to extremal SdS 
black hole is 1.5 (Ali and Rehman 2009). 

We consider the region — 5 < z < 5 for wave analysis assuming that 
event horizon is at z = 0. We take this region to explain waves near horizon 
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only for convenience. Since the flow variables exhibit large variations in 
the region — 1 < z < 1, we ignore it and solve dispersion relation for two 
meshes, i.e., —5 < z < — 1 and 1 < z < 5 (corresponding to near and far 
electromagnetic radiation zone). A complex dispersion relation ( Das 2004) 
is obtained by solving the determinant of the coefficients of constants of the 
corresponding equations of the rotating non- magnetized plasma. The real 
part of the determinant yields a quartic equation in k 

A^k 4 + A 2 {z, u)k 3 + A 3 {z, u)k 2 + A 4 {z, u)k + A 5 (z, u) = (3.30) 

which gives four real roots. A cubic equation in k is obtained from the 
imaginary part 

B^k 3 + B 2 (z, u)k 2 + B 3 (z, u)k + B 4 (z, u) = (3.31) 

which yields three real roots. The first and second root of the real part show 
wave propagation only in the region — 5 < z < — 1 while the third and fourth 
root exhibit waves in the region 1 < z < 5. The roots of the imaginary part 
indicate wave propagation in both meshes —5 < z < — 1 and 1 < z < 5, i.e., 
region near the event horizon and outer end of magnetosphere respectively. 

The wave vector, refractive index, its change with respect to angular 
frequency, group velocity and phase velocity lead to the wave properties of 
the SdS black hole and properties of Veselago medium. These are shown in 
Figures 1-10 by using real values of k in Eqs. fl3.30p and (13.311) . 

It is given that dispersion is normal if phase velocity is greater than the 
group velocity, otherwise anomalous (Achenbach 1973) or equivalently disper- 
sion is normal if change in refractive index is positive, anomalous otherwise. 
We see from figures that some waves move towards the event horizon and 
some move away from the horizon. The dispersion is normal in Figures 5-7 
and 10 while it is anomalous in the whole region of Figure 8. The following 
table classifies the regions of normal and anomalous dispersion in Figures 1-4 
and 9. 
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The results deduced from these figures can be expressed in the following 
table. 



Table I. Direction and refractive index of waves 



Fig. 


Direction of Waves 


Refractive Index (n) 






n 


< 1 and decreases in the region 


1 


Move towards the event horizon 




_5 < z < —lA, < oj < 2.6 
with the decrease in z 






n 


< 1 and increases in the region 


2 


Move away from the event horizon 




-2.94 < z < -1.1,0 < oj < 3.5 
with the decrease in z 






n 


< 1 and increases in the region 


3 


Move towards the event horizon 




1 < z < 1.9, 1.6 < co < 2.7 
with the decrease in z 






n 


< 1 and decreases in the region 


4 


Move outwards from the event horizon 




1 < z < 5, < oj < 3.8 
with the decrease in z 






n 


< 1 and decreases in the region 


5 


Move towards the event horizon 




-5 < z < -4, 5 < oj < 8 
with the decrease in z 






n 


< 1 and decreases in the region 


6 


Move away from the event horizon 




1.8 < z < 5,5.9 < oj < 10 
with the decrease in z 






n 


< 1 and decreases in the region 


7 


Move away from the event horizon 




-5 < z < -2.6,0 < oj < 4.7 
with the decrease in z 






n 


< 1 and decreases in the region 


8 


Move away from the event horizon 




1 < z < 5,6.3 < oj < 10 
with the decrease in z 






n 


< 1 and decreases in the region 


9 


Move towards the event horizon 




-5 < z < -3.4,8.1 < oj < 10 
with the decrease in z 






n 


< 1 and decreases in the region 


10 


Move toward the event horizon 




1.4 < z < 5,0.8 < oj < 5.4 
with the decrease in z 
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Surface Plot of Have Vector Surface Plot of refractive Index 




Surface Plot of D[n,w] Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 1: Dispersion is normal and anomalous in the region 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Figure 2: Normal as wells as anomalous dispersion occur at random points 
in the region. 
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Surface Plot of Wave Vector Surface Plot of refractive Index 




Figure 3: Normal and anomalous dispersion of waves is observed. 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Surface Plot of D[n,w] Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 4: Random points of normal and anomalous dispersion are found in 
the region. 
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Surface Plot of Wave Vector Surface Plot of refractive Index 




Surface Plot of D[n,w] Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 5: Whole region admits normal dispersion. 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Surface Plot of D[n,w] Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 6: Dispersion of waves is normal throughout the region. 



15 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Figure 7: Waves disperse normally in the whole region. 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Surface Plot of D[n,w] Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 8: Region shows anomalous dispersion. 



16 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Figure 9: Waves exhibit both normal and anomalous dispersion. 



Surface Plot of Wave Vector Surface Plot of refractive Index 




Surface Plot of D[n,£U] Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 10: Dispersion is found to be normal in the whole region. 
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Table II. Regions of dispersion 



Fig. 


Normal dispersion 


Anomalous dispersion 


1 


-5 < z < -1.3,0 < u < 10 
-1.2 < z < -1,0 < co < 1 


-1.28 <z< -1.2,3 < oj < 5 
-1.28 < z < -1.2,7 < uj < 10 


2 


-5 < 2 < -4, < cj < 10 
-4 < 2 < -3,0 < cj < 1 
-1.5 < z < -1,0 < u < 1 


-3 < z < -2.5,7 < w < 8.5 
-2.5 < z < -2,7 < w < 8.4 
-2< z< -1.5,2.5 < u < 3 


3 


3.9 < 2 < 5,0 < w < 10 


1.2 < z < 3.6,0 < w < 10 


4 


1.5 < z < 5,0 < u < 10 
1 < « < 1-5,0 < w < 2 


1 < z < 1.1,3 < u < 4 
1.1 < 2 < 1.3,5 < oj < 7 


9 


-5 < z < -1,0 < w < 4.5 
-4 < z < -3,4.7 < to < 5 


-5 < z < -4,8.2 < u < 10 
-4 < z < -3,8.5 < cj < 10 



4 Plasma Flow With Rotating Magnetized 
Background 

Here plasma is supposed to be rotating and magnetized. The magnetic field 
and velocity of fluid are assumed to lie in xz-plane. The corresponding 
perturbed Fourier analyzed GRMHD equations, i.e., Eqs. fl2.23p - (12.291) are 
given in Section 2. 

4.1 Numerical Solutions 

We take the same assumptions for the lapse function, velocity and specific 
enthalpy as in the previous section. Further, we assume ^ = 2 with u = V 

and V F = 1 in Eq.(j2HID so that A = 1 + ^±^. 

Here we also consider the region —5 < z < 5, 0<w<10 and investigate 
the wave properties in meshes —5 < z < — 1 and 1 < z < 5 From Eqs. 02.24"]) - 
f )2.25p . it follows that c 5 = 0. Consequently, we obtain dispersion relation 
whose real part is 

A^z)^ + A 2 (z, u)k 3 + A 3 (z, u)k 2 + A 4 (z, u)k + A 5 (z, u) = (4.1) 

giving four imaginary roots. The imaginary part of the dispersion relation 

B 1 {z)k? + B 2 (z, u)k A + B 3 (z, co)k 3 + B 4 (z, co)k 2 + B 5 (z, u)k 
+Bz(z,u) = (4.2) 
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yields five roots of k out of which one is real and four roots are complex. The 
real root indicates wave propagation in both meshes, i.e., —5 < z < — 1 and 
1 < z < 5 shown in Figures 11-12. This shows that waves move towards the 
event horizon. Also, it is obvious from figures that dispersion is normal as 
well as anomalous at random points. 

The following tables show the results obtained from these figures. 



Table III. Direction and refractive index of waves 



Fig. 


Direction of Waves 


Refractive Index (n) 


11 


Move towards the event horizon 


n < 1 and decreases in the region 
-5 < z < -2.1,0 < u < 10 
with the decrease in z 


12 


Move towards the event horizon 


n < 1 and increases in the region 
1 < z < 1.8, 1.5 < oj < 4 
with the decrease in z 



Table IV. Regions of dispersion 



Fig. 


Normal dispersion 


Anomalous dispersion 


11 


-5 < z < -4.9, 1 < u < 10 
-4 < z < -3,2.5 < u < 10 
-3 < z < -1,2.8 < u < 10 


-4 < z < -3.6, 1 < u < 1.5 
-3.7 < z < -3.5,1.9 < w < 2.1 
-3.5 < z < -3.35,1.8 < u < 2.2 


12 


1<2<2,2<u;<10 
2 < z < 4.5,5 < u < 10 
4.1 < z < 4.5,5.1 < u < 10 


1 < z < 2,0.8 < u < 1.1 
4 < z < 4.5,1.4 < co < 4 
4 < z < 4.5,1.4 < u < 4.5 



5 Summary 

This paper deals with the study of isothermal plasma wave properties in 
magnetosphere of SdS black hole in a Veselago medium. The ADM 3 + 1 
formalism has been used to formulate the GRMHD equations for this unusual 
medium. We have applied linear perturbations to the GRMHD equations and 
have obtained their component form with the assumption that plasma flows 
in two dimensions. Finally, we have obtained dispersion relations for the 
rotating (non- magnetized and magnetized) background. 

For the rotating non-magnetized background, waves move towards the 
event horizon shown in Figures 1, 3, 5, 9 and 10 while waves are directed 
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Surface Plot of wave Vector 



Surface Plot of Refractive Index 




Surface Plot of D[n,£j] Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 11: Normal and anomalous dispersion at random points. 



Surface Plot of wave Vector Surface Plot of Refractive Index 




Figure 12: Dispersion is normal as well as anomalous at random points. 
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away from the event horizon in Figures 2, 4, 6, 7 and 8. The dispersion 
is found to be normal as well as anomalous at random points in Figures 1, 
2, 3, 4 and 9. The Figures 5, 6, 7 and 10 show normal dispersion while 
8 admits anomalous dispersion in the whole region. The Figures 11 and 
12 indicate that waves are directed towards the event horizon for rotating 
magnetized plasma. It is clear from these figures that region admits normal 
and anomalous dispersion at random points. 

We know that the refractive index is always greater than one in the usual 
medium, while it is less than one for the Veselago medium. Here we have 
found that the refractive index is less than one and increases in small regions. 
The phase velocity is greater than group velocity for both non-magnetized 
and magnetized backgrounds. These are prominent aspects of the Veselago 
medium which confirms the presence of this unusual medium for both rotat- 
ing (non-magnetized and magnetized) plasma in SdS black hole. 

It is interesting to mention here that in a recent work (Sharif and Mukthar 
2011a, 2011b) for isothermal plasma on Schwarzschild black hole, there does 
not exist waves for the rotating magnetized plasma. However, we have seen 
wave propagation in SdS black hole for this case. Here waves admit normal 
dispersion at most of points while for the schwarzschild black hole, most of the 
waves disperse anomalously. Thus it can be concluded that more information 
can be extracted from magnetosphere by inclusion of the de-Sitter patch in 
the Schwarzschild spacetime. It would be interesting to extend this analysis 
for hot plasma which is in progress. 

Appendix 

The Maxwell equations, the 3 + 1 GRMHD equations for the SdS spacetime 
are given in this appendix. The Maxwell equations for such a medium are 



The GRMHD equations for the SdS spacetime in Rindler coordinates turn 




(Al) 
(A2) 



(A3) 



V x B = - 



(A4) 
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out to be (Sharif and Mukthar 2011a, 2011b) 

— = -V x («V x B), (A5) 

V.B = 0, (A6) 

do <9V 

+ (aV.V)po + Po7 2 V.^ + po7 2 V.(«V.V)V 

+ Po V.(aV) = 0, (A7) 
{(/W + ^)Sa + PotfVlVj - l-BiBM— + V.V)V* 

~ h BiB * )yi * yk + Po7V ^f + (V - VM 

= -p oAi7 2 a 4 -p n +i-(V x B)iV.(V x B) - J-^aB) 2 ,, 

B* - -L[B x {V x (V x («V x B))}],, (A8) 

+ V.V)W) - -% + W(V.a) + WIV.V) 

a at a at 

__(V x B).(V x - — (V x B).(B x -^-) 

+— (VxB).(Vx«B)=0. (A9) 
Aira 
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